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Abstract 

A model is discussed in which an electric field induces quantum nucle- 
ation of soliton-antisoliton pairs in a pinned charge or spin density wave. 
Coulomb blockade prevents pair creation until the electric field exceeds 
a sharp threshold value, which can be much smaller than the classical 
depinning field. We calculate the vacuum state energy and expectation 
value of the phase <f>, which is treated as a quantum scalar field. We find 
that the phase can also be much smaller, below threshold, than predicted 
by classical "sliding" density wave models. 
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Charge and spin density waves are examples of spontaneous symmetry break- 
ing, in which pairs of electrons and holes condense into a new ground state. 1 A 
charge density wave (CDW) forms in a linear chain compound when the elec- 
tronic charge density becomes modulated: p (x) = p (x)+p 1 cos [2kpx — (f> (x, t)], 
where 4> is the phase. A spin density wave (SDW) has a modulated spin den- 
sity, AS (x) = AS cos [2kpx — <j> (x, t)], and is equivalent to two out-of-phase 
CDWs for the spin-up and spin-down subbands. Although pinned by impuri- 
ties, a density wave (DW) can transport a current when an applied field exceeds 
a threshold value Ep- The simplest model of a pinned DW is a sine-Gordon 
(s-G) model, 2 in which the "vacuum states" are the 2im minima about which a 
DW would oscillate at the pinning frequency uq in the absence of dissipation. 
Finite wavelength modes, or phasons, can propagate at the phason velocity, 
Co = p~ 1 t 2 vp, where p = Mp/m e is the Frohlich mass ratio. 

The Hamiltonian for a pinned DW, including Coulomb interactions, can be 
written as: 



H 



= /^|^ + ^(^) 2 + ^§[l-cos^]+« B [^-fl] 2 |, (1) 



where D = — — (per spin per chain), 3 and the canonical momentum den- 

sity is given by 7r = Ddt4>. The last term in Eq. (1), introduced by Krive and 
Rozhavsky, 4 is the electrostatic energy density \e (E int ± E) 2 due to the applied 
field E and the internal fields E int ~ — (<^/27r) E* generated by phase variations. 
Here, E* represents the internal field produced by a 2ir soliton-antisoliton (S- 
S') pair (Fig. 1(a)), up = \sA (E*/2n) is the electrostatic energy density, 
9 = 2-kE/E* is proportional to the applied electric field, and <fr is measured 
with respect to its value at x = ±oo. The above Hamiltonian is equivalent to 
the bosonic form of the massive Schwinger model, as discussed by S. Coleman. 5 
A topological soliton (antisoliton) of width = c /ui n ~ lpm (>> A^) car- 
ries a charge =pe* = ^fn'e per spin-chain, where n' is the condensate fraction. 
A soliton-antisoliton (S-S') pair can nucleate by quantum tunneling when an 
electric field is applied. 6 In a related picture, a pair is created when a quantum 
soliton Zener tunnels through a "Bardeen pinning gap." 7 ' 3 

The electrostatic (Coulomb) energy of the S-S' pair introduces a sharply 
defined threshold field for pair creation. Figure 1(a) shows an S-S' pair, dis- 
playing the phase <fi an d excess charge density bp = — (e*/27r) d x <j> as func- 
tions of position. The pair, analogous to a parallel plate capacitor of sep- 
aration L and cross-sectional area A per chain (per spin), produces an in- 
ternal field of magnitude E* = e*/eA. The enormous dielectric constant, 
e ~ 10 8 £q, may include an intrinsic contribution Edw from the pinned DW, 
as well as a substantial contribution e s due to screening by the normal car- 
riers. When an external field E is applied, the difference between the elec- 
trostatic energy of a state with a pair and that of the vacuum is given by 



AU = -eAL 
2 



(E±E*f - E 2 



-E* ±E 
2 



Note that AU is positive 
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when \E\ < ^E* , so conservation of energy forbids the vacuum to produce a 

pair for fields less than a threshold value Et = ^E* (corresponding to 6 = ir). 

The threshold voltage across a region of length L will thus be e*/2C (where 
C = eA/L), 3 by analogy to Coulomb blockade in tunnel junctions. Above 
threshold, the S-S' nucleation and annihilation events will become correlated in 
time, by analogy to time-correlated single electron tunneling (Fig. 1(b)). 8 

This model has at least two important consequences. First, the threshold 
field for pair creation, Et (= e* /2eA), can be much smaller than the classical 
depinning field in the s-G model, E c = 2TrDLU$/e*, provided ue < Duuq. Sec- 
ond, the phase displacement <j> below threshold can also be much smaller than 
either the classical s-G prediction, <p = sin -1 (E/E c ), or the substantial phase 
displacements predicted by more realistic classical models. 9,10 A growing body 
of evidence demonstrates that the CDW phase displacements below threshold 
are extremely small for a wide range of temperatures in NbSe3 and TaS3— a 
fact which cannot readily be interpreted classically. This evidence includes: (1) 
the observed 2° displacement of CDW phase in NbSc 3 , as determined by NMR 
experiments, 11 when E = (3/4) Et, (2) the complete absence of any increase in 
CDW dielectric response ("critical polarization") below threshold in NbSc3 12 
and TaS3, 13 in serious contradiction with classical predictions, 9 ' 10 and (3) most 
recently, the observed absence of any change in CDW wavevector Q near the 
contacts below threshold. 14 

The primary goals of this paper are, treating the phase as a scalar quantum 
field, to calculate the vacuum state energy, e, and expectation value of the 
phase ((p) in the metastable regime. In order to evaluate e and ((f)), we follow 
Coleman 15 and employ a variational method, choosing the trial states to be 
coherent states. Such a choice is especially appropriate for a density wave, 
whose actual ground state (of the 3-D system) is, in fact, a coherent state. 
A CDW, for example, consists of a condensate of electron-hole pairs coupled 
to a condensate of 2fcp phonons, which are bosons. Such a ground state is 
characterized by a well-defined expectation value (</>}. It is pointed out in Ref. 
[15] that a piece of every term in the perturbation series expansion is effectively 
included by this variational approach. 

Rescaling (ct,x) = (x^x 1 ) — > — (x^x 1 ), the Hamiltonian Eq. (1) now 

Dc a 



reads: 



H 



with 



N„ 



00 Y 
5x° J 



dx 1 



VI 



dxH{x), (2) 



where 



U(4>) =lu 2 [1 -cos0] + -ml ((/> ■ 

,1/2 



of 



(3) 



-0^0' 



mo = [2-DcqUb] and N mo means normal order with 
respect to mo- The parameter m can be interpreted as a "free field" mass in 
the absence of the washboard pinning potential. 



3 



For the theory defined by Eq. (2), all ultraviolet divergences can be removed 
by normal-ordering in the interaction picture. However, variational calculations 
are carried out in the Schrodinger picture. Nevertheless, as pointed out in Ref. 
[15], one can still define a normal order N m with respect to an arbitrary mass 
parameter m. Following Ref. [15], the Hamiltonian density can be re- written 
as: 



H = N ma [Ho + U (</>)] = N m [Ho + U {<j>, m)} 



(4) 



where Ho is the kinetic piece of the Hamiltonian Eq. (2) and U (</>, m) is defined 
by: 



U (4>, m) = cxp 



8tt n \m 2 ) d4> 2 



U{<j>) + 



8tt 



(5) 



Now, it is straightforward to calculate the expectation value of the Hamiltonian 
density on the vacuum states appropriate to free fields of arbitrary mass m, 
which is given by the free coefficient of Eq. (5). Obviously, this class of states 
is not rich enough to be a good trial for the variational method. In order to 
ensure that {(f>) is well defined and to allow nonzero values of {(/>), one introduces 
a family of coherent states, |£ (x) ,p(x)), labeled by two arbitrary functions of 
space: £ (x) and p(x). Such coherent states have the property that, for any 
function of tt (x) and <f>(x), F (tt (x) ,(f>(x)), the expectation value is given by: 



N m F i 



(tt (x) ,j>(x))\ti(x),p(x))=F(p (x) ,£(x)). (6) 



These are used as trial states to minimize the expectation value (H) . Using Eq. 
(4) and the property Eq. (6) we have then: 

(£ (x) ,p(x) \H (x)\£ (x) ,p{x)) = \ V (xf + i (d^ (x) f + U (£ (x) , m) . (7) 

Choosing p = and £ (x) = £(= ((f))) = const, we minimize the first two terms. 
The energy density becomes: 



1-1 — 
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1 

8?7 
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8tt 



(8) 



Minimizing e with respect to m and £, one obtains the following self-consistency 
equations: 
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(9) 



61 = 0. 
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We have also derived the above self-consistency equations by extending the mean 
field approximation first developed by Glimm, Jaffe, and Spencer 16 for quartic 
interactions and then extended by Imbrie 17 to polynomial interactions. Details 
of our calculations will be reported elsewhere. 18 

Equations (9) enable one to determine both the mass-ratio m/mo and phase 
displacement £ = (</>) self-consistently in the metastable regime. Figure 2 shows 
the resulting plots of (cf>) vs. 9 (where 9 is the normalized applied electric field) 

for several different values of the parameter r = , which represents the 

DuJ o 

ratio of the electrostatic energy to the pinning energy The solid lines through 
the symbols represent the lowest energy states. Recall that S-S' pair creation 
can occur only when 9 > it (i.e. E > Et — E*/2). If 6 = n was taken to 
represent the classical depinning field (i.e. if Et was assumed to be equal to 
E c ), then the phase below threshold would be given by (<j>) = sin- 1 [E/E c ] = 
sin -1 [9 c /ir] according to the classical sine-Gordon model (as indicated by the 
x'ed line in Fig. 2). However, when r << 1, the phase displacements can be much 
smaller than these classical s-G predictions, right up to the critical value 9 = tt 
above which an instanton transition can occur without violating conservation 
of energy. (Note that Et << E c if r << 1.) The square in Fig. 2 represents 
the experimental result reported in Ref. [11], which is consistent with our 
calculations, provided we take the pinning energy Dujq to be much larger than 
the electrostatic energy ue (t « 0.02 in this case). When Dlo\ >> ue the 
threshold field for pair creation Et will be much smaller than the classical 
pinning field E c . The observed small phase displacements indicated by Refs. 
[11-14] arc all consistent with Et being small compared to E c . 

Figure 3 shows plots of energy density obtained using Eg. (8) vs. 9 for a 
case where the pinning energy dominates over the electrostatic energy (i.e. for 
r = 0.2). From Fig. 3, one can immediately see that the instanton transition is 
forbidden unless 9 > tt. When the applied electric field is below threshold (9 < 
7r), pair creation is prevented by Coulomb blockade. One can gain additional 
insight into the problem by considering the case when the applied field is equal 
to its critical value, i.e. 9 = tt. Figure 4 shows plots, obtained from Eqs. (9), 
of (</)) and the mass ratio m/mo a s functions of the parameter r for the case 
where 9 = tt. Note the bifurcation that occurs when r « 0.87. Thus, more than 
one metastable state will exist only when r is less than this critical value (i.e. 
when the pinning energy is greater than the electrostatic energy). 

The transition temperatures for DW formation can be quite high, well above 
200 K in a number of materials. The coupling between parallel DW chains and 
resulting 3-D coherence must therefore suppress thermal excitations of the S- 
S' pairs, whose energy would be extremely small if only a single transverse 
degree of freedom was considered. Ref. [3] provides heuristic scaling arguments 
on how thermal excitations might be suppressed without suppressing the S- 
S' tunneling amplitude for an intermediate range of coupling strengths, where 
the CDW electrons remain delocalized in real space. The situation here may be 
analogous to Josephson tunneling, where individual Cooper pairs tunnel through 
an insulating barrier even though the condensed pairs comprise a condensate 
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with a single thermal degree of freedom. A complete theoretical treatment of 
DW transport should consider many interacting scalar fields 4> n and, perhaps, 
generalize the concept of coherent Josephson-like tunneling. The implications of 
a collective quantum mechanism of DW depinning are potentially profound and 
far-reaching, so further work is warranted to experimentally test the predictions 
here, and to further develop a DW depinning model based on the principles of 
quantum field theory. 
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Fig.l (a) A materialized soliton-antisoliton pair, showing the position- de- 
pendent phase <j) (x), excess charge density 5p (x) = — (e* /2n) d<p/x, and inter- 
nal field E* — e*/eA. (b) The time-evolution of the phase (f)(x,t), illustrating 
the nuclcation and subsequent annihilation of three pairs during the first cycle. 
Nucleation of additional pairs is blocked in the lightly shaded regions (due to 
the internal field, -E*) until all pairs have annihilated, thus leading to time- 
correlated pair creation and annihilation. 

Fig.2 (</>) vs. 8 for r = 0.02 (stars), 0.2 (circles), 0.87 (dashed line), and 
1.5 (diamonds). The solid lines through the symbols indicate the lowest energy 
states. The x'ed line represents the classical sine-Gordon prediction (defining 
8 = 8 C = nE/E c ), while the square represents the experimental results obtained 
by Ross et al. [11] 

Fig. 3 Energy vs. electric field (8) for r = 0.2, showing the two branches 
corresponding to ((f)) « (stars) and (</>} w 2ir (circles). 

Fig. 4 ((f>) vs. r (left hand axis, circles), and mass ratio m/m vs. r (right- 
hand axis, stars) for 8 = n. Note that the system bifurcates when r « 0.87. 
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